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ABSTRACT 

We present observations of the CO[J=3^2] emission towards two massive and infrared luminous 
Lyman Break Galaxies aX z = 3.21 and z — 2.92, using the IRAM Plateau de Bure Interferometer, 
placing first constraints on the molecular gas masses (Mgas) of non-lensed LBGs. Their overall proper- 
ties are consistent with those of typical (Main-Sequence) galaxies at their redshifts, with specific star 
formation rates ~ 1.6 and ~ 2.2 Gyr~^, despite their large infrared luminosities (Lir ~ 2-3 x 10^^ 
Lq) derived from Herschel. With one plausible CO detection (spurious detection probability of 10"'^) 
and one upper limit, we investigate the evolution of the molecular gas-to-stellar mass ratio (Mgas/M») 
with redshift. Our data suggest that the steep evolution of Mgas/M* of normal galaxies up to z ^ 2 
is followed by a fiattening at higher redshifts, providing supporting evidence for the existence of a 
plateau in the evolution of the specific star formation rate aX z > 2.5. 
Subject headings: 



1. INTRODUCTION 

Until recently, molecular gas observations in the dis- 
tant Universe were restricted to very luminous, rare ob- 
jects, while the gas content of normally star forming 
galaxies - that fall within the so-called "Main Sequence" 
relation between their star formation rates and stellar 
masses (e.g., Brinchmann et al. 2004; Noeske et al. 2007; 
Elbaz et al. 2011; Magdis et al. 2010a) - at z > 2 was 
unknown. However, in the last few years the first molec- 
ular gas surveys of normal, massive-star-forming galaxies 
at z = 1 - 2 were performed (Daddi et al. 2008; 2010; 
Tacconi et al. 2010), revealing that such typical distant 
star-forming galaxies were gas rich. Lyman Break Galax- 
ies (LBGs) are the most common population of normal 
star- forming galaxies at z ^ 3 (e.g., Steidel et al. 1999, 
Adelberger et al. 2000) and are ideal targets to extend 
this study to even higher redshifts. 

Multi- wavelength studies of LBGs have provided ex- 
tensive information on various physical properties of 
these objects such as stellar masses and dust attenua- 
tion (e.g., Papovich et al. 2001, Rigopoulou et al. 2010, 
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Magdis et al. 2008, 2010a, b). However, their gaseous 
component still remains largely unconstrained, as cur- 
rently the only two LBGs with molecular gas measure- 
ments are two gravitationally lensed systems: MS1512- 
cB58 (Yee et al. 1996) at z = 2.72 and the LBG 
J213512. 73-010143 -"the Cosmic Eye"- (Small et al. 
2007), at z = 3.07. Baker et al. (2004), Coppin et al. 
(2007), and more recently Riechers et al. (2010) have 
reported CO detections for cB58 and Cosmic Eye, pro- 
viding evidence of the existence of a sizeable cold gas 
reservoir in LBGs. The fact that these two galaxies 
are intrinsically faint and probably not representative of 
massive, z 3 galaxies, coupled with the uncertainties 
introduced by lensing, implies that more data and di- 
rect CO observations of massive, normal high— z galaxies 
are indispensable for an improved understanding of these 
systems. 

In this letter, we present CO measurements for two 
z ~ 3, non-lensed LBGs in the northern field of the Great 
Observatories Origins Deep survey (GOODS-N), selected 
to have properties representative of normal galaxies at 
this redshift. We couple these observations with data 
from the Herschel Space Observatory (Pilbratt et al. 
2010), obtained as part of the GOODS- Herschel program 
(PI D. Elbaz), to study in detail their infrared proper- 
ties and place constraints on their gas content and on the 
evolution of the gas to stellar mass (Afgas/M,) evolution 
from z = 2 to z = 3. We assume 51^ — 0.3, Ho = 71 km 
sec-i Mpc^i, Q.X = 0.7 and a Chabrier IMF. 

2. SAMPLE AND OBSERVATIONS 

The two galaxies of this study (M23 and M18) were 
originally optically selected ([/„, G, i? with R < 25.5) by 
Steidel et al. (2003) in the GOODS-N field. Rest-frame 
ultraviolet (UV) spectroscopy - which was also used to as- 
certain the absence of AGN signatures (i.e. strong high- 
ionization emission lines) - provided redshifts, Zspec = 
3.216 and z^pec = 2.929, respectively, determined from 
Ly-a (among other) absorption lines. The sources 
benefit from extensive multi-wavelength (UV to radio) 
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Fig. 1.— ACS V-band (5" x 5"), MIPS 24 /im (40" x 40"), 
PACS 160 Mm (40" x 40") and SPIRE 250 /im (200" x 200") cut- 
out images of M23 (top) amd M18 (bottom). The cyan crosses are 
centred at the IRAC 3.6 ^m position of the sources. 



coverage, including ground based Un,G, R, H, J, K ob- 
servations and photometry from the Advance Camera 
for Surveys (ACS, B,V,i,z), Infrared Array Camera 
(IRAC) and the Multi-band Imaging Photometer (MIPS 
24 /im) on board Spitzer and the Very Large Array 
(1.4 GHz). The GOODS-77ersc/ieZ program provides 100- 
and 160 fim imaging with the Photodetector Array Cam- 
era and Spectrometer (PACS, Poglitsch et al. 2010) and 
250- , 350- and 500 /xm with the Spectral and Photomet- 
ric Imaging Receiver (SPIRE, Griffin et al. 2010 ). Her- 
schel photometry was performed using the source extrac- 
tion point-spread function (PSF) fitting code galfit (Peng 
et al. 2002), guided by 24 /xm priors (for more details see 
Daddi et al. in prep). Both sources have a > Scr detec- 
tion in at least one Herschel band. Images of the two 
sources in several bands, are shown in Fig. 1. 

We observed the CO[J=3^2] (t/^cst = 345.796GHz) 
transition towards the M23 and M18, using the Plateau 
de Bure Interferometer (PdBI). At z = 3.216 and z = 
2.929, this line is redshifted to 82.0199 and 88.0112 GHz 
respectively. Observations were carried out under good 
3 mm weather conditions in D array configuration with 
5 antennas during May- June 2011, resulting in 3.8h and 
3. Oh on source time after flagging of bad visibilities. Data 
reduction was performed using CLIC in GILDAS. The 
noise level is 0.34 mJy/beam over 300 km s~^ in both 
datasets, and the FWHM of the circularized synthetized 
beam is about 6" for both. 

3. RESULTS 

We searched for CO emission at the position of the 
LBG targets (M23, RA: 189.2612457, DEC: 62.2405968 
and M18, RA: 189.1837769, DEC: 62.2196693) and de- 
tected possible emission lines at S/N = 4.0 and 3.3 for 
M23 and M18, respectively. The first is exactly at phase 
center while the latter reaches 3.6cr if we allow for a small 
offset of 1.4" (only 20% of the synthesized beam). In both 
cases this peak emission is within the velocity range al- 
lowed by the la optical redshift uncertainties of 600 km 
s~i and 1000 km s"^ respectively, for M23 and M18. This 
is much smaller than the >12000 km s~^ velocity range 
sampled by the 3.6 GHz WIDEX correlator. In order to 
quantitatively assess the likelihood of chance/spurious 
detection of similar signals we extracted spectra from 
our datasets at 10000 fixed positions chosen within twice 
the primary beam but excluding the central LBG region. 
For each of this spectra we consider 600 km s^-^ (1000 
km s~^) subsets at a time for M23 (M18) and system- 
atically look for maximally positive signals within this 



range. From this we estimate that the chance probabil- 
ity to have S/N > 4.0 for M23 is 1.4 xlO'^. Based on i) 
the fairly small chance of the signal being spurious, ii) the 
fact that the observed signal is spatially coincident with 
the position of the source (Fig. 2 left) and iii) that the 
emission is centred at the expected frequency within la 
of the optical spectroscopic uncertainty (Fig. 2 middle 
and right), we consider M23 to be a plausible detection, 
and derive a CO[3-2] luminosity of (1.60 ± 0.4) x lO^'' K 
km s~^ pc^. For M18, we only place a 4cr upper limit of 
< 2.2 X 10^° K km s^^ pc^ (assuming a FWHM of 400 
km s~^), as the detected signal is weak, with a spurious 
probability of 3 x 10~^. 

4. DISCUSSION 

4.1. Physical Properties of the Galaxies 

We combine all available photometric data, ranging 
from rest-frame UV to observed radio 1.4GHz, to con- 
struct the full SEDs of the sources (Fig. 3). We fit 
the infrared part of the SEDs using the Chary & Elbaz 
(2001) templates as well as the Main-Sequence and star- 
burst SEDs, described in Magdis et al. (2011b), to de- 
rive estimates for the infrared luminosity of the sources. 
All templates result in similar estimates, indicating that 
both sources have ULIRG-like luminosities. In particular 
we find iiR = (2.9 ± 1.6) x 10^^ fo^. ^23 and Lir = 
(3.0 ± 1.1) x 10^2 L© for M18. Given the low S/N of 
the infrared data, we also employ the radio detections 
(SiA = 20.1 ± 5.2/iJy for M23 and 15.3 ± 4.8 /xJy for 
M18, Morrison et al. 2010) for an independent estimate 
of the LiR of the systems. In particular we convert the 
radio fluxes to Li.4GHz, assuming and slope oi a = —0.8 
and based on the Condon (1992) IR-radio correlation, we 
flnd Lm = 4.8 x lO^^ L© for M23 and Lir = 1.9 x lO^^ 
Lq for M18, in reasonable agreement with those derived 
based on SED fltting. 

Fitting the rest-frame UV-to-near-IR part of the spec- 
trum with the Bruzual & Chariot (2003) model SEDs 
assuming a constant star formation history, yields a stel- 
lar mass of 2.0 x 10" M© for M23 and 9.5 x 10^° Mq for 
M18. Based on the best flt models the sources are found 
to be moderately obscured with similar extinction values 
( E{B — V) ~0.25 ). Furthermore, using the extinction- 
corrected UV-luminosity L1500 and adopting the Kenni- 
cutt (1998) relation we derive dust corrected star forma- 
tion rates that correspond to Ljr « 1.6 x 10^^ Lq and 
1.3 X 10^2 foj. ]y[23 and M18 respectively. These UV 
based Lir estimates are ~2 times lower compared to the 
"true" LiR derived from the IR data. Nevertheless, this 
discrepancy is considerably lower than what is found for 
local ULIRGs and high— z starbursts, and closer to that 
found for normal galaxies (e.g., Goldader et al. 2002, 
Elbaz et al. 2007, Magdis et al. 2010, Rigopoulou et al. 
2010). 

Given that the far-IR, radio and UV Ljr estimates 
agree within the uncertainties, we choose to average the 
three indicators and subsequently derive the best SFR 
measurements. With SFR = 310 Mg yr'^ and 210 Mq 
yr~^ for M23 and M18 respectively our analysis yields 
sSFR ^ 1.6 Gyr~^ and ~ 2.2 Gyr"-'^ for the two sources. 
These values are very close to the characteristic, mass de- 
pendent sSFRms - derived based on the SFR- A/* Main 
Sequence of z 3 LBGs presented by Magdis et al. 
(2010) - with sSFR/sSFRms '-I -2 and 0.9 for M23 and 
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Fig. 2. — CO observations of M23. left; Contours (from ± 2a and in steps of Icr) of CO[3-2] emission from M23 overlaid on IRAC 
3.6 /im images (40" size), middle: Spectra of CO [3-2] emission binned in steps of 75 km s~^ for the full range of velocity covered by our 
observations, right: Zoomed version of the middle panel, presenting only the 25% of the total available spectral range. The red color 
indicates the regions where positive emission is detected. These regions have been used to derive total integrated fluxes. The vertical 
dashed linos indicate the velocity range that corresponds to the redshift uncertainty (Icr) as derived by optical spectroscopy (see Figure 3). 




Fig. 3. — Observed rest-frame UV to radio SEDs of M23 (left) and M18 (right). For the SEDs we present the detections (purple squares) 
along with 5(t upper limits (purple arrows). The rest frame UV-near-IR portion of the data is overlaid with the best fit BC03 model (black 
solid line), while the mid-IR to radio is shown with the best fit CEOl model (orange solid line) and the best fit models of SB (green) and 
MS (black) templates from Magdis et al. (2011). Inset panels, Keck/LRIS rest-frame UV spectra of the sources (courtesy of A. Shaplcy). 
Red lines correspond to the best fit LEG spectral templates of Shapley et al. (2003). 



M18 respectively. Consequently, despite the high SFRs 
of our galaxies that are much higher than those of most 
common galaxies at z ~ 3, their overall properties 
are consistent with those of normal rather than to those 
of starburst galaxies that tend to have substantially ele- 
vated SFRs for their stellar masses. 

4.2. Molecular Gas Mass 

CO measurements can be used to derive estimates of 
the molecular gas mass of a galaxy (e.g., Solomon & Van- 
den Bout 2005). However, the conversion factor from CO 
luminosities to the molecular gas mass, Mgas — aco x 
L'qq, introduces a large uncertainty in this derivation. In 
particular Downes & Solomon (1998) showed that aco is 
a factor of ~ 6 smaller for local ULIRGs, than for local 
spiral galaxies with average values of 0.8 (but also see Pa- 
padopoulos et al. 2012) and 4.6, respectively. Quite nat- 
urally, it has been a common practice to apply the same 
dco value for both local and high— z ULIRGs. Neverthe- 
less, recently, there has been substantial evidence (both 
theoretical and observational), that the aco value for 



high— z Main Sequence galaxies, is comparable to that 
of local spirals and that it depends on the metallicity 
rather than the infrared luminosity of the source (e.g., 
Magdis et al. 2011, Genzel et al. 2011, Narayanan et al. 
2011, Daddi et al. 2010). In light of these results, and 
given that our LBGs are Main Sequence galaxies, we will 
attempt to make an educated guess of the aco value for 
our LBGs, based on the aco ~ metallicity relation for 
normal galaxies (Genzel et al. 2011). 

Since we lack direct measurements of the metallicity 
for the two sources, we employ the fundamental metal- 
licity relation of Mannucci et al. (2010) that relates 
the SFR and the stellar mass to metallicity. We de- 
rive a nearly solar metallicity for both sources, (i.e., 
12+log{0/H) = 8.7) that based on the Genzel et al. 
(2011) relation would give aco^^ ^ 3.6. A similar 
value would be derived based on the local aco - Z re- 
lation of Leroy et al. (2011). We also note that this 

The units of acd Mqpc^'^ (K km s^-"^)"^, are omitted from 
the text for brevity. The quoted values include helium. 
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TABLE 1 
Observed and Derived Properties 



Source 2° mI Lir" /co[3-2] S/N det LcxP /gas 

M0 L0 Jy km K km s^^ pc^ 

M23 3.216 (2.0±0.3) x lO" (3.1±1.1) x 10^^ 0.36±0.08 O (3.20±0.8) X 10^" (1.15±0.3) xlO" 0.36 

M18 2.929 (9.5±1.4) x 10^° (2.1±1.3) x 10^^ <0.52 - <4.40 x lo'-" <1.58 xio'^'^ <0.62 

Notes: 

a: A(z) = 0.004 and 0.006 for M23 and M18 respectively 
b: averaged Lir from UV, far-IR and radio estimators 
c: from CO[3-2] emission line, derived after applying rsi — 0.5 
d: assuming aco— 3.6, including He 



value is very close to the one derived by Magdis et al. 
(2011) and Daddi et al. (2010), for normal BzK galax- 
ies, using dust mass measurements and dynamical con- 
straints, respectively. We convert luminosities derived 
from C0[J=3— >2] to L^q (defined relative to the funda- 
mental C0[J=1— »^0] transition), adopting an excitation 
correction of rsi = 0.5 based on the CO spectral line en- 
ergy distribution measured by Dannerbauer et al. (2009) 
for a BzK galaxy at z ~ 1.5. Then based on the above we 
derive M^^ = (1.15 ± 0.30) x 10" x (aco/3.6) Mq for 
M23 and an upper limit of < 1.58 x 10" x (aco/3.6) 
M0 for M18. 

Focusing on M23, we infer a gas consumption timescale 
(Tgas = -Wgas/SFR) of ^ 0.4 Gyr, very similar to that of 
normal disks at z ^ 1-2 by Daddi et al. (2010) and Gen- 
zel et al. (2010) and a factor of ~ 7, higher than that of 
local ULIRGs. This suggests that it can maintain its star 
formation activity for much longer time scales, than those 
expected for rapid, merger-driven bursts, highlighting a 
fundamental difference in the star-formation mode of lo- 
cal ULIRGs and high— z Main Sequence galaxies with 
comparable infrared luminosities. We note that although 
this result is subject to several uncertainties, it is unlikely 
to be driven by systcmatics in the adopted aco value. 
By simply examining at direct observables, L'f~,Q and Lm , 
M23 has a star formation efficiency (SFE = Lir/Lco) ^ 
95 L© (K km s~^ pc^)~^, similar to that of high— z disks 
(Daddi et al. 2010). This is considerably lower when 
compared to the two Icnscd LBGs in the literature for 
which CO observations arc available, the Cosmic Eye and 
the cB58 (Riechers et al. 2010), with Lir/L^q~260 and 
~600 Lq (K km s~^ pc^)^^ respectively. Furthermore, 
the two lenscd LBGs, exhibit an excess in their sSFR by 
a factor of ~6 with respect to Main Sequence galaxies at 
this redshift (Riechers et al. 2010). The elevated sSFR 
and the high SFE of the lensed LBGs, typical character- 
istics of merger driven star-bursting systems (e.g., Daddi 
et al. 2010, Elbaz et al. 2011), arc in striking contrast 
with the "normal" star forming mode of M23. Evidently 
a larger sample of non-lensed 2; ~ 3 LBGs is essential to 
drive this investigation. 

4.3. The evolution of the Gas-to- Stellar Mass Ratio 

Currently, one of the most notable debates between 
theoretical predictions and observational evidence, is the 

evolution of sSFR beyond z ~ 2. In particular, sev- 
eral observational studies (but not all, see de Barros et 
al. 2012 for example) indicate that the sSFR evolves 
steeply up to z ~2, and then remains roughly constant 
throughout the redshift range z = 2 — 7 (e.g., Gonzalez 
et al. 2010; Daddi et al. 2010). On the other hand, cur- 
rent theoretical approaches fail to reproduce this plateau, 
predicting a gradual increase of sSFR at higher redshifts 



(e.g., Dekel et al. 2010, Weinmann et al. 2011, Dave, et 
al. 2012). Given that the gas mass of normal galaxies is 
intimately linked with their star formation with Mgas « 
SFRO-81 g^i^ redshifts (e.g., Daddi et al. 2010), the evo- 
lution of gas-to-stellar mass ratio (Mgas/M*) should be 
closely linked to the evolution of the sSFR, albeit with 
a shallower slope. Indeed, since sSFR=SFR/M*, it im- 
plies that at a fixed stellar mass, Mgas/A^* oc sSFR"-^-^. 
Therefore, Mgas estimates can be used to trace the evo- 
lution of sSFR in an independent way. 

While restricted to only two galaxies, we can place con- 
straints on the evolution of Mgas/M,^, or equally of the 
molecular gas fraction (/gas = A'/gas/[Mg„.,+M*]), of nor- 
mal galaxies up to z ^ 3 by combining our data with re- 
sults from previous studies of normal galaxies at various 
redshifts. In Figure 4, we present estimates of Mgas/M*, 
for normal galaxies in the local universe from Leroy et al. 
(2008), at z - 0.4 from Geach et al. (2011), at z - 0.5 
from Salmi et al. (2012) (in prep.), and at z ^ 1 — 2 from 
Daddi et al. (2010) and Tacconi et al. (2010), including 
only sources with > 10^° Mq. For our LBGs we find 
that Mgas/M* w 0.6 (M23) and Mg^/M^ < 1.6 (M18). 

To test the two competing scenarios regarding the evo- 
lution of sSFR at z > 2.5, we plot in Figure 4 the evo- 
lution of Mgas/A'/* for the case of a gradually increasing 
sSFR as (l-|-z)^-^ at all redshifts, driven by the dark mat- 
ter specific accretion rate as predicted by the theoretical 
work of Dave ct al. (2012). We also plot the observation- 
ally motivated tracks of Sargent et al. (2012b in prep), 
that are built assuming Mgas oc SFR°-^^ at all redshifts 
and an increase of sSFR as (1 -|- z)^'® up to z ~ 2.5 fol- 
lowed by a flattening towards higher redshifts, as well 
as the evolutionary tracks of Reddy et al. (2012) for L^, 
galaxies at z ^ 2 — 7. Our Mgas/M* estimate and upper 
limit, appear to be consistent with a flattening of the 
evolution, providing supporting evidence for a plateau in 
the evolution of sSFR too, at z > 2.5. We note that an 
anti-correlation between Mgas/M* and M^, is expected 
given that (i) the sSFR is not mass-invariant but de- 
clines with M* and (ii) the relation between Mgas and 
SFR is not linear. Therefore, the fact the our LBGs are 
2 to 4 times more massive than the comparison sample 
(that has (M*) « 5 x IO^^Mq) could introduce a bias in 
our result. However, rescaling the Mgas/M* of all galax- 
ies in Figure 4 to M* = 5.0 x 10^° Mq by assuming 
M^^/M^ oc M-O "^ (e.g., Daddi et al. 2010, Saintonge et 
al. 2011, Magdis et al. 2012, submitted), does not affect 
our conclusion. Finally, the flattening of the evolution 
would be even more pronounced if the observed signal 
from both of our sources were treated only as upper lim- 
its and is further observationaUy supported by the non 
detection of CO emission from z ~ 5 LBGs (Davies et 
al. 2010). 
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Fig. 4. — Evolution of A/gas/M, and /gas with redshift. Literature data are: orange circles from Leroy et al. (2008), empty black circles 
from Geach et al. (2011), blue squares from Daddi et al. (2011) and Salmi et al. (2012) in prep., green squares from Tacconi et al. (2010). 
The 2 ~ 3 LBGs of this study are shown with filled red squares. M18 is presented as a 4(t upper limit. Dashed lines depict the tracks of 
Sargent et al. (2012b in prep), for the case of Af* = 5.0 X 10^° Mq, M* = 1.0 X lO" Mq and M, = 2.0 X 10^^ Mq. Cyan lines show the 
evolutionary tracks of Reddy et al. (2012), for galaxies with a final stellar mass M* = 5.0 X IQl" Mq at z = 2.3 for the case of a constant 
sSFR = 2.4 Gyr-l at z > 2 and for the case of sSFR = 2.4 Gyr'^ up to 2 ~ 3.0 and sSFR = 5.0 Gyr'^ at 2 ~ 4. The solid orange line 
depicts the predicted evolution of Mgas/M« for the case of an increasing sSFR oc (1 + 2)2-9 at all redshifts (Dave et al. 2012). 
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